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Abstract 



We investigate the applicability of the moduli space approximation in the- 
ories with unbroken non-Abelian gauge symmetries. Such theories have mass- 
less magnetic monopoles that are manifested at the classical level as clouds 
of non-Abelian field surrounding one or more massive monopoles. Using an 
SO (5) example with one massive and one massless monopole, we compare the 
predictions of the moduli space approximation with the results of a numerical 
solution of the full field equations. We find that the two diverge when the 
cloud velocity becomes of order unity. After this time the cloud profile ap- 
proximates a spherical wavefront moving at the speed of light. In the region 
well behind this wavefront the moduli space approximation continues to give 
a good approximation to the fields. We therefore expect it to provide a good 
description of the motion of the massive monopoles and of the transfer of 
energy between the massive and massless monopoles. 
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I. INTRODUCTION 



The moduli space approximation (MSA) [jTl is a useful tool for studying the low-energy 
dynamics of Bogomolny-Prasad-Sommerfield (BPS) magnetic monopoles and other sys- 
tems that allow families of static multisoliton solutions. In this approximation the full field 
dynamics is replaced by that of the small number of collective coordinates that span the 
moduli space of static solutions. This neglects possible distortions of the inner structure 
of the solitons, and thus effectively treats them as point particles. The approximation is 
expected to be valid when the soliton velocity f <^ 1, with corrections suppressed by powers 
of V. This suggests that it might break down in certain gauge theories where, as we describe 
below, some BPS monopoles are effectively transformed into massless particles Indeed, a 
study by two of us |^ of monopole scattering in one such theory encountered anomalies that 
suggest a breakdown of the MSA. In this paper we describe some analytic and numerical 
calculations that shed light on this issue. 

The MSA can be motivated by recalling some features of the spectrum of small oscilla- 
tions about a multisoliton solution. First, there are a number of zero- frequency modes that 
are handled by the introduction of the collective coordinates. Next, for each elementary field 
of mass iTLi in the theory there is a continuum spectrum beginning with minimum energy 
E = rrii. Finally, there may be some discrete nonzero eigenvalues. The last two components 
of the spectrum correspond to scattering states and bound states, respectively, of the mas- 
sive quanta in the presence of the solitons. If the soliton kinetic energies are much smaller 
than the m^, excitation of these modes should be suppressed. In the MSA one assumes that 
this suppression is complete and takes the collective coordinates to be the only dynamical 
degrees of freedom. If, as is the case for BPS solutions, the static energy is independent of 
the collective coordinates, the effective Lagrangian for the collective coordinates has only 
kinetic energy terms. These define a metric on the moduli space, and the classical evolution 
of the system is equivalent to geodesic motion with respect to this metric. 

The energetic arguments underlying the MSA apply fairly straightforwardly to theories 
with only massive fields. However, when massless fields are present, as is the case in theories 
with magnetic monopoles, radiation of low-energy quanta is always energetically possible. To 
justify the MSA, one must show that this radiation is suppressed at low monopole velocities. 
The simplest case to consider is that of an SU(2) theory broken to U(l), with gauge coupling 
e, a massive vector boson of mass m, and a monopole mass M ~ Airm/e'^. The problem of 
radiation in this theory was first addressed by Manton and Samols p. For well-separated 
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monopoles the amount of U(l) "electromagnetic" radiation^ can be estimated by treating 
the monopoles as point magnetic charges moving along the trajectories specified by the 
MSA. Standard electromagnetic techniques show that the total dipole radiation is of order 
Mf ^, while higher multipoles are suppressed by additional powers of v. (For the case of two 
monopoles, the dipole moment vanishes identically, so the quadrupole term dominates and 
the total radiation is of order Mv^ .) These arguments break down if the monopoles approach 
closely enough that their cores actually collide and overlap. However, since the core radius 
is of order m~^, the dominant effect should only be on modes of wavelength ^ m~^, with 
quanta of energy ^ m, and these can neglected because of energetic considerations. These 
somewhat heuristic arguments have been placed on a more rigorous footing by Stuart 0. 

Now consider BPS monopoles in a theory with a gauge group G of rank k > 1. For generic 
vacuum expectation values of the adjoint Higgs field the symmetry is broken to the maximum 
Abelian subgroup U(l)'^'. There are then k distinct topological charges. Associated with each 
is a fundamental monopole p[, carrying one unit of that charge, that can be realized as an 
embedding of the SU(2) unit monopole in the SU(2) subgroup defined by one of the simple 
roots of G. All higher- charged BPS solutions can be understood as multimonopole solutions 
containing appropriate numbers of the various species of fundamental monopoles. As long 
as the fundamental monopoles all remain massive, the arguments used to justify the MSA 
in the SU(2) theory can be readily extended to this more general case. 

The situation is more complex if the Higgs field is such that the unbroken gauge group 
has a no n- Abelian component. When this happens, the masses of some of the fundamental 
monopoles vanish. These massless monopoles, which can be viewed as the counterparts of 
the massless non-Abelian gauge bosons, cannot be realized as isolated classical solutions. 
Instead, they are manifested in classical multimonopole solutions as "clouds" surrounding 
one or more massive monopoles |Q. Both Abelian and non-Abelian magnetic fields are 
present within the cloud region, but outside the cloud the non-Abelian components are 
suppressed by powers of r relative to the Abelian Coulomb magnetic fields. For static 
solutions the energy is independent of the size of the cloud. 

One might well expect to encounter some difficulties with the MSA in the case. First, 
the presence of massless non-Abelian gauge fields gives more possibilities for low-energy 
radiation. Second, the charged fields that are the source for this radiation are not confined 
to a core region of size ~ m~^, as in the case of Abelian breaking, but are instead present 
throughout the cloud region, which can be arbitrarily large. Finally, the vanishing mass of 
some of the monopoles seems inconsistent with the low-velocity regime that is required for 
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the MSA. 



This last point needs some elaboration. Because the massless monopoles do not have well- 
defined positions, their velocities are also not well-defined. However, there is an associated 
quantity with dimensions of velocity, namely the rate of expansion of the cloud region. This 
need not be unity, as an ordinary massless-particle velocity would be. On the other hand, 
even if this "cloud velocity" is initially small, there is no guarantee that it will remain small, 
or even bounded. For example, in Ref. 0, the MSA was used to study solutions containing 
two massive monopoles and a massless cloud in an SU(A^) theory. It was found that the 
massive monopoles tend toward a constant velocity at large times, but that the cloud size 
varies roughly as kt'^, where k is proportional to the kinetic energy associated with the cloud. 
Thus, no matter how small this energy, the MSA predicts that the cloud eventually expands 
faster than the speed of light. 

The prediction of such superluminal expansion in a relativistic field theory suggests a 
breakdown of the MSA. In principle, one could test this by solving the full time- dependent 
field equations and comparing with the MSA predictions. For arbitrary configurations of the 
SU(A^) monopoles this would be a fairly difficult undertaking. We therefore focus instead 
on another example that shares many of the essential features of the SU(A^) solutions but 
that is analytically more tractable. 

The example we consider occurs in a theory with an SO (5) gauge symmetry broken to 
SU(2)xU(l) by an adjoint Higgs field. There are two species of fundamental monopoles, one 
massive and one massless. The static BPS solutions containing one of each are all spherically 
symmetric and contain a massive core of fixed radius surrounded by a cloud of arbitrary 
radius a. Suppressing the dependence on the other collective coordinates, we denote the 
corresponding fields by A^^^{r, a) and $^^^(r, a). 

If we work in the center-of-mass frame and restrict ourselves to solutions with vanishing 
SU(2) and U(l) electric charges, the cloud radius is the only time-dependent collective co- 
ordinate. Solving the Euler-Lagrange equation following from the moduli space Lagrangian 
gives a solution amsaI^)- The MSA consists of assuming that the time-dependence of the 
fields arises solely from that of a, and that the latter is given by QmsaI^)- Explicitly, 



To test this approximation, we start with a configuration that is of this form at t = 0; 
i.e., that satisfies 



<'^(r,t)=<^^(r,aMSA(t)), 




(1.1) 



A^(r,0) = A 



(r,ao) 
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$(r,0) = $^PS(r,ao), 

MfS(r,ao) 



da 



<l>(r,0) = do (1-2) 

with dots denoting time derivatives. Since we expect the MSA to be most rehable when 
velocities are small, we take the initial cloud velocity dg -C 1. We also assume that the 
initial cloud size is much greater than the radius of the massive core. 

These initial conditions determine a solution of the full field equations that can be com- 
pared with the MSA solution of Eq. ( p..l| ). Although the spherical symmetry implied by the 
initial conditions simplifies the field equations somewhat, one is still left with a large number 
of coupled differential equations. However, the fact that the departures from the MSA are 
expected to occur primarily outside the massive core leads to a further simplification that 
reduces the problem to a single second-order equation. Analysis of this equation leads to a 
semi-quantitative understanding of the breakdown of the MSA. To test these ideas in detail, 
however, numerical methods are required. 

In Sec. II we review the SO (5) theory in detail and describe the moduli space of BPS 
solutions. In Sec. Ill we set up the initial value problem that we will be considering. Here we 
analyze the field equations and obtain predictions for the expected departures from the MSA 
solution. In Sec. IV we describe the numerical solution of the field equations, and compare 
it with our expectations. In Sec. V we review our results and discuss their implications for 
the use of the MSA. 



II. SO(5) MONOPOLES 

The generators of SO (5) may be chosen to be two commuting operators Hi and H2 that 
generate the Cartan subalgebra, together with eight raising and lowering operators -Ei/j, 
where the i^j are the roots shown in Fig. |l|. (We fix the normalization of the gauge coupling 
by taking the long roots to be of unit length.) The vacuum expectation value of the adjoint 
Higgs field can always be chosen to lie in the Cartan subalgebra, and hence to be of the 
form $00 = hiHi + /i2-f^2 = h ■ H. Furthermore, h can always be chosen so that h ■ (3 and 
h ■ 7 are both positive. The two species of fundamental monopoles are then associated with 
the roots f3 and 7 and have masses 

Mf3 = ^ = — h-/3, Mj = J: = _h-7. 2.1 

1^ e f3 e ' e 7^ e 



4 



The corresponding electrically-charged elementary gauge bosons have masses 

m^ = eh-f3, = eh ■ . (2-2) 

If h ■ /3 and h ■ 7 are both nonzero, the unbroken gauge group is U(l)xU(l). If instead 
h is orthogonal to 7, the symmetry is only broken to SU(2) xU(l), and the 7-monopole 
mass, like m^y, vanishes; this is the case we wish to consider. For convenience, we will define 
V = h - a = 2h - (3. 

at 




Y 



FIG. 1. The root diagram of S0(5) 



Index theorems show that there is an eight-dimensional moduli space of BPS so- 
lutions containing one massive /3- and one massless 7-monopole. Three of the collective 
coordinates can be chosen to be the center-of-mass coordinates X, while four others can be 
taken to be the angles a, P, 7, and x that specify the global SU(2)xU(l) orientation of 
the solution. Variations of these seven variables are equivalent to the actions of symmetry 
transformations. The dependence of the solutions on the last collective coordinate, which 
we take to be the cloud parameter a, is less trivial. 

To display this dependence, we need some notation. We assemble the generators into a 
triplet 



E-a), 



V2 



{Ecx - E_cx) 



U a) 



a H: 



(2.3) 



a second triplet ta(7), defined by analogous expressions involving 7; and a quartet assembled 
into a 2 X 2 matrix 
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M = V2i(^^ ^ ^ I . (2.4) 

The spacetime fields can then be written as 

A, = Af ^ ■t{a) + Af ■ t(7) + tr Af M , 

$ = $(1) • t(a) + $(2) . t(^) + tr ^^^^M . (2.5) 

With the SU(2) x U(l) symmetry breaking that we are assuming, the asymptotic Higgs field 
lies in the subgroup generated by the t{ct). The field components labeled by superscripts 
1,2, and 3 transform as three singlets, a triplet, and a complex doublet, respectively, under 
the unbroken SU(2). 

In Ref. it was shown that there is a static solution of the form 

^"(1) = ea^mrmA{r) , $^1) = faH{r) , 

A^(3) = r,F(r) , $(3) = -tIFir) . (2.6) 



Here 



and 



A{r) = (2.7) 

smh evr er 



Hir) = V coth evr (2-8) 

er 



are the functions appearing in the SU(2) unit BPS monopole, while 



F(r) = L(r, af'^ (2.9) 

^ ' ^cosh(ez;r/2) ^ ' ' ^ ' 



and 

G{r) = A{r)L{r,a) , (2.10) 

with 

L(r, a) = [1 + (r/a) coth(ewr/2)]"^ . (2.11) 

At the center of this solution is a massive core of radius ~ 1/ev. Surrounding this is a 
non-Abelian cloud region whose size is set by the parameter a, which can take on any positive 
real value. Inside the cloud, in the region 1/ef <^ r <^ a, the non-Abelian vector potential 
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A falls as 1/r and yields a Coulomb non-Abelian magnetic field. In the region r ^ a, the 
non-Abelian magnetic charge is canceled by the cloud, and A-^^ falls more rapidly, as 

Given the explicit form of the solution, it is a straightforward matter to obtain the metric 
on the moduli space, and thus the moduli space Lagrangian M 



>MS — —MryX? -\- , 

^2 Mae^' 



27r 



a2 



X H — j{ ha a +sin a (3 +(7 + cosa/5)' 



(2.12) 



III. FIELD EQUATIONS AND ANALYTICAL PREDICTIONS 

Our goal in this paper is to understand how well the solutions given by the MSA approx- 
imate the solutions of the full field equations. To simplify matters, we will restrict ourselves 
to solutions with vanishing linear momentum and SU(2) xU(l) charges. The Euler-Lagrange 
equations that follow from Eq. ( p.l2|) then imply that the only time-dependent collective 



coordinate is a, given by 

ausKit) = kit + /\f (3.1) 
where k and A are constants. The associated kinetic energy is 

E = 2f^ = ^. (3.2) 

The MSA approximation to the fields is obtained by substituting Eq. ( |3.1| ) into Eq. ( |1 . 1| ) 
We want to compare this with the solution of the full field equations. Because of the spherical 
symmetry of the static BPS solutions, our initial conditions imply that the solutions will be 
spherically symmetric for all t. Even with this taken into account, the most general Ansatz 
for the fields would involve a large number of functions of r and t. However, we expect 
that the deviations from the MSA arise primarily from the dynamics of the massless fields, 
rather than from the massive fields that are confined to the monopole core. Furthermore, we 
expect these deviations to manifest themselves primarily at large distances. Let us therefore 
choose some radius f that is large compared to the monopole core radius, but much smaller 
than the initial cloud size; i.e., 1/ev <^f <^ oq. We will assume that the MSA gives a good 
approximation to the fields at r <^ f for all t. (We will see that our numerical results are 
consistent with this assumption.) We then only need to solve the field equations for r > f, 
subject to the boundary condition that the solution match onto the MSA solution at r = f. 

This leads to considerable simplification. The massive SU(2) doublet fields A^f) and ^'^^^ 
fall exponentially outside the monopole core and can therefore, to a good approximation, 
be set equal to zero in the region r > f. Once this is done, the singlet and triplet fields are 
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decoupled from each other. Since the former are independent of a, we can concentrate on 
the latter. 

Requiring spherical symmetry and positive parity implies that the triplet fields must be 
of the form 

Af^=fJ{r,t), 

$-(2) =f„[G(r,t)-ir(r,t)]. (3.3) 

J can be set equal to zero by a time-dependent gauge-transformation. The Euler-Lagrange 
equations of the theory then reduce to 

~G + G" + -C -^G- —G(2G + K)- 2e\G + K)G^ = (3.4) 



-k + K" + -K' - —K - —K{2K + G)- 2e\G + K)K^ = (3.5) 



/y ijn^ ijn 

where dots and primes denote differentiation with respect to t and r, respectively, and 
exponentially small quantities have been ignored. 

The initial conditions of Eq. (|1.2| ) imply that both K and K vanish for all r at t = 0. 
It then follows from Eq. ( p.5|) that K will continue to vanish for all time. (One can also 
verify that small deviations from if = do not produce growing perturbations.) Defining 
g{r,t) = —erG{r,t), we can then write 

--g + g"~^gil-g)' = 0. (3.6) 

This has the static solution 

gMs{r,a) = —^ (3.7) 
r + a 

which corresponds to the r ^ 1/ev limit of the BPS solution Eq. ( p.lOj ). Note that g ^ 1 
far inside the cloud, while far from the cloud g tends to zero. 
The initial conditions of Eq. ( |1.2|) reduce to 

g{r,t = 0) = fi'Ms('^,ao) , 



g{r,t = 0) = ao 



da 

In addition, we have the boundary condition 



a=ao 



(r + aof 



(3.8) 



9{r, t) = gusir, aMSA(t)) = _ ^msa(^) _ ^^ g-^ 

r + aMSA{t) 
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The time dependence of this boundary condition imphes that the total energy in the region 
r > f will not be conserved. Physically, this is to be expected. Although the total energy of 
the static solutions is independent of a, the division of this energy between the core and the 
cloud varies with cloud size. Hence, even within the MSA there will be energy flow across 
the surface r = f. 

We will find it necessary to use numerical methods to solve Eq. ( ^.61 ). We will describe 
these in the next section. Before doing so, let us consider what type of departures from the 
MSA we might expect. We can distinguish two different ways in which the MSA might be 
expected to fail. First, the solution may depart from the moduli space by the excitation of 
low-frequency modes; i.e., by radiation. In the present context, this can occur because the 
time-dependent fields implied by the MSA act as a source for non-Abelian radiation. To 
make this more precise, let us write g{r, t) = gusif^, dusAit)) + 5g{r, t). If 5g <^ gus, we can 
linearize Eq. (|3.6|) to obtain 



2 

+ V — ^(1 - fl'Ms)(i - ^9Ms)Sg = gus 



{r + omsa)^ {r + omsa)^ 
e^Er{r - Somsa) 



(3.10) 



47r {r + a^sAf 
where E is the cloud kinetic energy given by Eq. ( |3.2| ) 

The right hand side of this equation is concentrated in the cloud region, r < a. (For this 
initial analysis of the problem we ignore the 1/r tail.) Hence, we can view gu_s as the source 
for a radiation field 5g that propagates outward from the cloud at the speed of light. In 
principle, we could solve for 5g in terms of this source by using Green's function methods. 
However, because of the presence of the terms involving (^ms; we have not been able to obtain 
a closed form expression for the Green's function. 

Nevertheless, we can still estimate the rate of radiation; i.e., the power P that passes 
through a sphere of some fixed radius r ^ oq. From Eq. ( p.lO|) we see that 5G = 6g/{er) 
is proportional to eE. Since the radiation energy density is quadratic in derivatives of 6G, 
the power must be of the form 

Pr.e'E'f{d)r^E-f{a), (3.11) 
a 

with dimensional arguments showing that / cannot depend on a. It is also easy to see that 
/ cannot depend on e. Although further analysis would be required to determine the form 
of /(d) at small d, it is presumably of order unity when a is of order unity. 

The MSA must certainly break down by the time that the integrated power becomes 
comparable to the total energy of the cloud. Comparing Eqs. ( p.2|) and ( p.ll|) and noting 
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that the bulk of the energy loss occurs at large d, where / is of order unity, we see that this 
happens when t ~ a/d^ ~ k~^. This is also the time that Eq. ( |3.1| ) predicts a cloud velocity 
comparable to the speed of light, a second indication of the breakdown of the MSA. 

Because of its finite speed of propagation, the radiation described above can have no 
effect on the behavior of the fields at large distances; i.e., at r > t. However, there is a 
second type of breakdown of the MSA that occurs in this region. The MSA assumes that 
the system can be well described by a small number of time-dependent collective coordinates. 
This assumption can only be valid if retardation effects can be ignored, so that the system 
effectively reacts as a single unit to changes in the collective coordinates. This should be 
true within the cloud region as long as the characteristic time a /a is much greater than a; 
i.e., as long as the cloud velocity is much less than the speed of light. However, no matter 
how small d may be, we should expect departures from the MSA at distances r ^ a/a. 

We can examine this in more detail by using a large-distance expansion to solve Eq. ( p. 61) . 

Let 

= i:^- (3.12) 

n=l ' 

Substituting this expansion into Eq. ( |3.(j| ) gives a series of equations for the Cn(t). The first 
two of these require that Ci = C2 = 0. Integrating these with the initial conditions implied 
by Eq. ( |3.8| ) gives 

Ci = ao -\- dot , 

C2 = — Oq — 2aodot . (3.13) 



Hence, 



g{r, t) ^ gusir, ao + dot) . (3.14) 



In other words, the large-distance fields are close to the moduli space of BPS solutions, but 
with a cloud size that grows at the initial velocity do rather than at the accelerating velocity 
predicted by the MSA. 



IV. NUMERICAL RESULTS 

We now turn to our numerical results. As explained in the previous section, we assume 
that the MSA is valid within a region of radius f ^ 1/ef, and so only need to obtain a 
numerical solution for the region r > f. The field equations then reduce to the single equation 
Eq. ( p.6|) , which is to be solved subject to the initial conditions of Eq. (|3.8|) and the boundary 
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condition of Eq. ( |3.9| ). To fix the initial conditions we must specify the initial cloud size, oq, 
and the initial cloud velocity, do- Because the cloud size sets the characteristic length scale 
for the spatial variations of the field, discretization errors will be reduced if ^ r; for our 
solutions we take oq = 50f. (The actual numerical value has no direct physical meaning, 
since our approximations effectively set 1/ef, the only other physical length scale, equal to 
zero.) In order that we initially be in a regime where the MSA can be expected to be valid, 
do should be nonrelativistic; we choose do = 0.01. The parameters in Eq. are then 

= 2.5 X 10-^ao ^ and A = 200ao. 
We expect the MSA to break down at roughly the time that dMSA becomes of order 
unity; i.e., at t ~ tcrit, where 

tcrit = ^ = 2 X lO^ao . (4.1) 

According to Eq. (|0|), amsa would then be ~ 1/(4A;) = lO^ao- Hence, to see the breakdown 
of the MSA we must follow the evolution of the system until the cloud has grown by several 
orders of magnitude. In order to make the computational load more manageable, we proceed 
as follows. Starting at t = 0, we numerically evolve the system over the range f < r < rmax, 
where rmax = lOOOao. At r = f, we require that g match onto the MSA solution. To fix 
the behavior at the other end of the interval, we use an expansion such as that given in 
Eq. ( |3.12| ) to obtain an analytic approximation to g for rmax < r < oo. We then require that 
the numerical solution match onto this analytic approximation at r = rmax- As time goes 
on and the cloud expands, the spatial range of the numerical integration must be enlarged. 
To do this, we define a position ri/2{t) by the condition g{ri/2,t) = 0.5. When ri/2 grows 
to twice its original size, we replace rmax by r'^s.x = Sr^ax and at the same time double the 
step sizes in our integrations. The new initial data is given for f < r < rmax by the previous 
numerical results and for rmax < r < r'^^^ by the large- distance analytic approximation. We 
continue in this fashion, doubling the interval and step sizes as required. 

In Fig. we show our results for g{r) at six different times. For comparison, in each 
case we also show the predictions of the MSA. (Note that the distance scale changes from 
one panel to the next.) For the two earliest times, t = lOlao and t = IOSOoq, the MSA is 
clearly a good approximation to the exact solution. The only noticeable departure from the 
MSA appears to be due to the retardation effects at large distance; indeed, note the sharp 
bend in the data at r = t that is visible at t = IOSOoq and all later times. The third time 
shown is approximately tcrit/2. At this point the departures from the MSA are beginning to 
be significant, although the MSA still gives a fairly good qualitative picture. The times in 
the last three panels are all much greater than tcrit- In these we see the numerical solution 
departing from the MSA and evolving toward a step-like profile in which g falls from 1 to 
in a short interval near r = t. This step is essentially the radiation transformed into a 
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0.5 1 1.5 2 0.5 1 1.5 2 
r/t r/t 

(e) (f) 

FIG. 2. Plots of the numerical solution (solid line) and MSA approximation (dotted line) for 

g{r,t) at various times. Note that the distance scale varies from panel to panel. The panels 

correspond to (a) t = lOlao, (b) t = lOSOao, (c) t = 1.11 x lO'^ao, (d) t = 1.49 x lO^ao, (e) 
t = 5.45 X 10^ ao, and (f) t = 1.07 x lO^ao. 
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sharp wavefront moving outward at the speed of hght; it arises because the cloud that is 
the source of the radiation is now itself expanding at the speed of light. Although the step 
appears to be becoming narrower with increasing t, this is due to the change in distance 
scale from one panel to the next. Indeed, we see from Eq. ( |3.6| ) that far from the origin g is 
a function only of r — t, so that both the shape and the width of the step are constant at 
large time. 

Despite this clear breakdown of the MSA at large distances, we find that it remains a 
good approximation to the data at shorter distances. To illustrate this, in Fig. |] we plot 
the fractional deviation of g from the MSA prediction as a function of r/t for the times 
corresponding to the last three panels of Fig. For any fixed r/t < 1 this fractional 
deviation decreases with time. 



-0.2 



-0.4 - 



S -0.6 



-0.8 




-1 I ' ' ' ^ 

0.5 1 1.5 

r/t 

FIG. 3. Fractional deviation of g{r) from the MSA prediction. From inner to outer, the curves 
correspond to t = 1.49 x lO^ao, 5.45 x lO^ao, and 1.07 x lO^ao- Note that the distance scale 
increases with time. 

Another way to measure the deviation from the MSA is to define an r-dependent cloud 
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parameter a{r,t) by substituting the actual g{r,t) into Eq. ( |3.7| ). In other words, 

rg{r,t) 



a{r, t) 



l-gir,t) ■ 



(4.2) 



If the MSA were exact, a(r, t) would be independent of r and equal to aMSAit). In Fig. || we 
plot a(r, t) for several times ranging from AOa^ to SOao- At all these times a(r, t) is relatively 
flat near the origin, gradually falls, and then becomes approximately constant for r > t. 
The increase in the actual value in the latter region is very close to being linear in time, as 
predicted by Eq. (|3]T|). 




200 



FIG. 4. Plot of a{r,t). Reading from the bottom, the curves correspond to t = 40ao, 50ao, 
60oo, 70oo, and 80ao. The distances are given in units of oq. 

We show the corresponding data for some larger times in Fig. |^, although with a rescaling 
of axes to allow easier comparisons between the data at the four different times. Again we 
see that a{r,t) is roughly constant at smaller r, varying by no more than 10% for r < t/4. 
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FIG. 5. Plot of a(r, t). Reading from the top, the curves correspond tot = lOlao, 1.03 x 10^ oq, 

1.1 X 10^ ao, and 1.49 x 10^ oq. Note that the distance scale increases with time. 

V. CONCLUDING REMARKS 

In this paper we have investigated the rehabihty of the MSA when apphed to BPS 
monopoles in a theory with an unbroken non-Abehan symmetry. Whenever there are mass- 
less particles, the existence of non-BPS states arbitrarily close in energy to the BPS solutions 
raises the possibility that, no matter how small the kinetic energy, a time-dependent solution 
will deviate from the moduli space of static solutions to a sufficient extent to invalidate the 
MSA. For the case of Abelian massless fields, previous investigations have shown that the 
radiation processes that might give rise to such departures fall as a power of velocity and 
so can be neglected to leading approximation at low energy. We have found that this is not 
the case when the massless fields are non- Abelian. 

The time-dependent fields in the cloud region surrounding the massive monopoles act as 
sources for non-Abelian radiation. As in the Abelian case, radiation is suppressed at low 
velocities; in particular, the power radiated is proportional to the square of the cloud kinetic 
energy. The crucial difference lies in the duration of the radiation. The Abelian radiation 
in a multimonopole system falls off as the massive monopoles recede from each other. By 
contrast, the non-Abelian radiation in our system continues unabated as the cloud expands. 
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The net effect is that the field profiles depart sharply from the MSA prediction by the time 
that the cloud velocity approaches unity. 

Despite this, we do not find a complete breakdown of the MSA. It remains quite reliable 
as long as the cloud velocity is small. Even when the cloud has become relativistic, the 
MSA gives a good description of the fields well inside the cloud (say for r ^ t/4). If 
the massive monopoles within the cloud remain nonrelativistic, they will always be in the 
region where the MSA is valid. Hence, one should be able to use the MSA to describe both 
the motion of the massive monopoles and the transfer of energy between the massive and 
massless monopoles. Furthermore, the cloud does survive as a clearly identifiable object in 
the classical solution, although the smoothly varying profile of the static BPS solution is 
transformed into a relativistic wavefront. In a sense, the modifications to the MSA that we 
have found are the minimal ones consistent with the relativistic bound on velocity. 

Although we have focussed on the case where there is a strictly massless monopole, we 
should make a few remarks on the situation when one of monopoles is massive but much 
lighter than all the others. For the analogue of our SO (5) example Eq. ( |3.2| ) is replaced by 



where r is the separation of the two monopoles, /i is the reduced mass of the monopole pair, 
and tq = An/e^fi gives the size of the lighter monopole core when the monopoles are well 
separated. 

When r ^ ro, the two monopoles are overlapping, with the core of the lighter monopole 
enclosing the heavier monopole. This core is then much smaller than its natural size and 
has a radius approximately equal to r. In this regime r grows linearly with time, as in the 



massless case, and only tends to a constant value ^j2E/fi after the two cores separate. If 
E <^ fi, the two monopoles are always nonrelativistic, and the MSA is reliable. If instead 
E ^ fi, a. breakdown of the MSA, similar to that which we have seen in this paper, occurs 
while the two monopoles are still overlapping. We would expect this to lead to a deformation 
of the core profile of the lighter monopole so that initially, like the cloud in the massless 
case, it resembles a spherical wavefront expanding at the speed of light. However, when the 
core radius become of order ro, this wavefront presumably splits into a shell of radiation 
that continues to expand and a massive core that is left behind. The heavier monopole will 
eventually separate from this core, but not from the outgoing shell of radiation. 

For this two-body system the validity of the MSA is assured if the initial relative velocity 
is small, since the final configuration is nonrelativistic if the initial one is. This need not be 
true in a three-body system. An example occurs in an SU(3) theory broken to U(l)xU(l) 
with one of the fundamental monopoles being much lighter than the other. Consider a 
configuration containing one of the light and two of the heavy monopoles. Because there 




(5.1) 
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can be a net transfer of energy from the heavy to the hght monopoles, it is possible for the 
hght monopole to initially have f -C 1 but to emerge with a relativistic final velocity. As a 
result, the MSA fails at large times, as was first pointed out by Irwin |Tl[ . 



Finally, let us consider what happens as we move toward the strong coupling (large e) 
limit, where the duality conjecture leads us to expect the solitons and elementary particles 
to exchange roles. Even without taking into account the quantum corrections to the classical 
field profiles, we can see how the classical soliton picture is lost. For the massive monopoles 
this is signalled by the fact that the core radius ~ l/(e^M) becomes less than the Compton 
wavelength, so that a nonrelativistic monopole cannot be sufficiently localized for its solitonic 
properties to be directly observable. For the massless monopoles, the transition can be seen 
by focussing on the two phases in the evolution of the massless cloud. For a time interval 
of order tcrit ~ ^^e cloud behaves like a slowly expanding soliton, as predicted by 

the MSA. Beyond this time, it resembles a wave front expanding at the speed of light. In 
the strong coupling limit, tcrit-E' is much less than unity. The uncertainty principle then 
implies that the period when the massless monopoles behave like nonrelativistic solitons is 
unobservably short; instead, they will always appear as a wavefront of radiation moving 
at the speed of light. In retrospect, this makes it seem quite natural that, even in the 
weak-coupling regime, the cloud cannot be kept nonrelativistic. 
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